Parameter estimation for mixed states from a single copy 
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Given a single copy of a mixed state of the form p — Xpi + (1 — A)p2, what is the optimal 
measurement to estimate the parameter A, if pi and p2 are known? We present a general strategy to 
obtain the optimal measurements employing a Bayesian estimator. The measurements are chosen to 
minimize the deviation between the estimated- and the true value of A. We explicitly determine the 
optimal measurements for a general two-dimensional system and for important higher dimensional 



I. INTRODUCTION 

The estimation of quantum states is one of the basic 
primitives in quantum theory. This general task appears 
in several modifications and situations. When many 
copies of the state are available, one may try to obtain 
maximal information on the quantum state by employ- 
ing state tomography [l| , although the effort may become 
overwhelming, if the dimension of the system increases. 
Furthermore, tomography is clearly not a viable way, if 
only few, or even only a single copy of the state are avail- 
able. 

Nevertheless, if some a priori information about the 
state is given, there are several estimation problems 
which can be meaningfully posed and solved even on the 
single copy level. The most prominent is state discrimi- 
nation, where one knows that the system is in one of sev- 
eral given states, and one needs to decide in which one. 
This problem has been discussed in several variations, ei- 
ther one aims for an unambiguous discrimination, or for 
a discrimination with a minimal error probability P, Q • 
Moreover, programmable state discriminators have been 
proposed [3|. Other problems considered are state esti- 
mation from several copies with collective or separable 
measurements [3| , the estimation of certain state param- 
eters (like the time in an unitary evolution) , as well as 
estimation of the state after a generalized measurement 
when the premeasurement state is unknown [6|. 

In this paper, we consider the following problem: Let 
us assume that an apparatus is given, which produces for 
a given input A the state 



p(A) = Api + (1 - \)p2 



(1) 



where pi and p2 are known. Let us further assume that A 
is determined by some well characterized random number 
generator, but we don't know its actual value. The task 
is now to estimate A from a single copy of p(A). 

Such estimation problems can occur in several realistic 
situations. For instance, one may consider some process 
of decoherence, where the single copy of p(A) describes 
the state of a single atom coupled to its environment, 



and the task is to estimate the rate of decoherence. In 
this case pi represents an initial state, while p2 is the 
state of thermal equilibrium, which the system assumes 
eventually. Another example may be that p(A) is the re- 
duced state of some multipartite pure state |-0(A)), where 
one tries to estimate A. We will discuss the decoherence 
example later in more detail. 

For the estimation one has to perform some measure- 
ment. The must general measurement is described by 
a positive operator valued measure (POVM) and we ask 
for the POVM that minimizes the expected deviation be- 
tween the true value of A and the estimated one. Here 
we take the viewpoint that the value of A is not exactly 
known beforehand but we may have information about it 
in form of a prior probability distribution. After a mea- 
surement this prior probability distribution can be up- 
dated according to Bayes theorem. Minimizing the devi- 
ation of the estimate from A expected from the posterior 
probability distribution leads to the so-called Bayesian 
estimator (cp. for example Q)- 

Bayesian estimators have been successfully applied in 
state discrimination problems. They are also appropri- 
ate in situations where a sequence of consecutive mea- 
surements is carried out and the estimate has to be up- 
dated after each measurement. This is for example the 
case in real-time monitoring of the dynamics of quan- 
tum systems. Using a Bayesian estimator it is possible 
to monitor oscillations of single qubits in real-time with 
high accuracy by means of a sequence of weak measure- 
ments d, Q . The Bayesian estimator is also intimately 
related to the continuous estimation of the wave function 
of a system with arbitrary-dimensional Hilbert space by 
means of continuous measurement Ahhough the 

search for the best measurement in connection with an 
Bayesian estimator is performed in this article for the 
particular task to estimate the parameter A in Eqn. ([1]), 
it may be carried out in an analogous way in different 
set-ups. 

For the case when p is a two-level system, i.e., a qubit, 
we determine the optimal measurement strategy for arbi- 
trary pi and p2. Interestingly, for the case of completely 
unknown A, it turns out that if pi and p2 do not have 



2 



the same purity, the optimal measurement does not com- 
mute with pi — P2, and hence does not coincide with the 
intuitive guess which one would make from the Bloch 
sphere picture. Finally, we discuss how the results from 
the qubit case can readily be used to solve this problem 
for important higher dimensional cases. 

As already mentioned, similar parameter estimation 
problems have in generality been studied in Ref. @ and 
general conditions on the information which can be ob- 
tained about the parameter have been formulated. Our 
aim, however, is to explicitly construct the optimal mea- 
surements. The knowledge of these explicit measure- 
ments may further be used to improve schemes for the 
observation of oscillations with small disturbance [1, Q . 

Our paper is organized as follows: In Section II we pose 
the problem in mathematical terms and derive a general 
condition on the optimal POVM. In Section III we show 
how this general condition can already reduce the set of 
POVMs which have to be considered. In Section IV we 
solve the problem for qubits. We first derive the optimal 
measurement, if pi and p2 are pure, then we consider 
the general case. We also discuss how one can apply 
the results for estimating decohercnce rates. Finally, in 
Section V we derive some results for higher dimensional 
systems. 



II. BAYESIAN ESTIMATOR AND CONDITION 
FOR AN OPTIMAL MEASUREMENT 

In this paper, we consider the following task. Let us 
assume that we have a single quantum system with d- 
dimensional Hilbert space TL. Its state is known to be a 
mixture of two states pi and /92, i-e.. 



px = Xpi + (1 - A)p2 



(2) 



where for the parameter A only some probability distri- 
bution is known. What is the measurement which leads 
to an optimal estimation of the unknown parameter A? 

The estimation of A has to proceed in two steps: 
First, a measurement is carried out on the system. The 
statistics of a general measurement can be characterized 
by a POVM and throughout this paper we denote by 
-p = {E^} a POVM with the effects Em- In a second 
step the parameter A is estimated by a number (the 
estimate) which depends on the result of the measure- 
ment TO. 

Let us first discuss the construction of the optimal es- 
timate gm for a given outcome m. A criterion for the 
optimality of the estimate g„i can be formulated as fol- 
lows: gm is optimal iff the cost function 

c(A,g„i) = (A - g™)^ (3) 

is expected to assume its minimum: 

£m ((A - /(A - gmfp{\\m)d\ = min . (4) 



Here the expectation value £ is taken with respect to the 
posterior probability distribution of A, which includes the 
information contained in the occurrence of the measure- 
ment result m: 



p(A|to) 



p(to|A)p(A) 
p(to) 



(5) 



where p{m) = J p{'m\X)p{X)dX is the probability to ob- 
tain the measurement result m averaged over all possible 
occurring states. 

Taking into account the linearity of the expectation 
value, one can directly verify that the optimal estimator 
gm is equal to the expected value of A: 



£m ((A - gm?) - £rn ((A)^) 



+ 9m : 



(fm(A) -5™)' + Var„,(A), (6) 



where 



Var„,(A) (A^ A)2 



(7) 



represents the variance of A. The right-hand side of ^ 
assumes a minimum for gm = £m (A). Such a value gm is 
also called a Bayesian estimate. 

Having derived the optimal estimate g,„ for a certain 
outcome to we can consider the optimal choice of the 
POVM. Now, the optimal measurement is the one which 
leads to the smallest expected costs averaged over all out- 
comes to: 



p{m)£m ((A - 5m)^) 



(8) 



Hence the effects {Em} have to minimize the mean vari- 
ance X]m?'('^)V^''ni('^)- Note that this variance, just as 
the expectation of the costs, is defined with respect to 
the posterior probability density p{X\m). 

For the sake of simplicity we choose in the following 
an equally distributed prior density p{X) = 1. We will 
see later in an example, that this imposes not a big re- 
striction. Using p(to|A) = tT[EmPx\ and the probability 
to obtain measurement result to 

p(to) = / p{m\X)p{X)dX 



tr 



tY[Em {Xpi + (1 - X)p2)]dX 



Em^{pi + P2) 



(9) 



we obtain for the posterior probability density in Eq. ([S]) 

tY[EmPx] 



p(A|to) 



tr [Em\{pl+P2)\ 



(10) 



In order to minimize the mean variance of A we have 
to first compute the expectation value of A and A^ with 
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respect to the posterior probability density: 

/q \tr[E„^px]dX 



Xp{X\m)dX - r . . 

tr[Em2(Pl+P2)\ 



_ tr [E^ + IP2)] 
~ 2tr[£;„i(pi +P2)] 

and in a similar calculation 



(11) 



/o 3tr [£;„i(/9i + P2)] 

The mean variance thus amounts to (see Eqns. ([71 [H])): 

™ 34^ tr[£;™i(pi+p2)] 

(12) 

We can thus summarize the main result of this section: 
Proposition 1. The optimal POVM V = {E^} is the 
one which maximizes 



7n L 



E^^ (pi + P2) 



(13) 



Note that using the fact that |pi + |p2 = 5(^1+^2) 
i(pi — P2) the quantity Q(7') may be rewritten as 

\2 



(tr[g,„(pi-p2)]) ^ 
3tr [i;^ (pi + P2) 



(14) 



which is manifestly invariant under the permutation of 
pi and p2- Further note that Eq. (fTTjl allows us to calcu- 
late the Bayesian estimate gm of A, given a measurement 
outcome associated with the effect operator Em-, since 

5m = i?m(A). 

The condition for an optimal measurement in case of 
a general prior probability distribution p{X) can be de- 
termined from the corresponding mean variance, which 
is of similar form as (fT^ . It can be derived as in Eqs. ([U 
[TUl [TTl [T^ with a general p(A) and reads: 



^p(m)Var„(A) 



= A2 



tr[£;„,(Api + (l-A)p2)] 



l-f)p. 



(15) 



where A" is the nth moment of p{X), i.e., A" = 
/q A"p(A)dA. Again, this may be rewritten in a way simi- 
lar to Eq. (|14p which is invariant under permutation of pi 
and p2 ■ Although we consider in the following an equidis- 
tributed A, thus p(A) = 1, we illustrate by an example 
in section IIVI that the results we obtain can readily be 
transcribed for the general case. 



III. PROPERTIES OF OPTIMAL 
MEASUREMENTS 

Now we show two facts about the POVM which mini- 
mizes the expression of Eq. (|12p : We first show that this 
POVM can be chosen to have pure effects, i.e. effects 
which are of rank one. Then, we show that this POVM 
must be a so-called extremal POVM. 

To show that the POVM which minimizes the expres- 
sion of Eq. (fT2)) has pure effects Ei , assume that we have 
a POVM where one effect, say Ei, is not pure. Then we 
can decompose Ei as Ei = E^ + Ef where E^ and E^ 
are linearly independent positive operators, which can 
serve as new effects. The new POVM with the effects 
E^ , El , i?2 , £"3 , . . . gives the same or a smaller value for 
the cost function. Indeed, from Eq. it follows that 
it suffices to show 



7B\2 



(16) 



where 



tr[E-{yi 



tr[£;^i(pi 



ip2)] and N° = Lr[£.^- ^ 
P2)\- Eq. p6)) can be straightforwardly verified, using the 
facts that Zi = Z^ + and Ni ^ + Nf^ . Thus, for 
our search for optimal measurement strategies, it suffices 
to consider POVMs with pure effects. 

There are further constraints on the optimal POVM 
V. These follow from the fact that the set of all POVMs 
is a convex set. Given two POVMs V'-^^ = {E^^ and 
75(2) — {Em^} with K outcomes, the convex combination 
V = pV^^^ + (1 -p)p(2) = {pEi^^ + {l-p)Ein^ is again 
a POVM with K outcomes. On the other hand there are 
POVMs which can not be expressed as a convex combi- 
nation of two different POVMs, these are called extremal 
[ll[. For our purpose, it is important to note that for 
convex combinations 



Qiv)<pQir^'y) + ii-p)Q{r^^^) 



(17) 



holds. This follows because it is true for each of the sum- 

mands of Q. This can be verified e.g. for the summand 

with the effect Ei = pE[^^ -I- (1 — p)e[^^ by replacing 

the effects E^ and E^ in inequality by pE[^^ and 
(2) 

(1 — p)El , respectively. Since a generic convex combi- 
nation on the right hand side of ([T7)) assumes only values 
less than Max{Q{r^^^), QiV^^^)}, only extremal POVMs 
can maximize Q and thus minimize the expected costs. 
We can summarize: 

Proposition 2. When maximizing Q in Eq. (fT^)) it suf- 
fices to consider extremal POVMs with pure effects. 

Indeed, we can require both conditions at the same 
time, since POVMs with pure effects, which are not ex- 
tremal, can be written as a convex combination of ex- 
tremal POVMs with pure effects. The two conditions on 
the POVMs are a priori independent. However, for spe- 
cial cases, certain relations are known. For instance, it 
has been shown that all extremal POVMs for qubits have 
pure effects and maximally four outcomes [TT| . This fact 
we will exploit in the following. 
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IV. OPTIMAL POVM FOR QUBITS 

111 this section we consider the problem for qubits, that 
is, pi and p2 are qubit states. The representation of these 
states on the Bloch sphere will enable us to determine the 
optimal POVM. We will first derive some general formu- 
lation of the optimization problem in terms of Bloch vec- 
tors. Then we will solve the problem for the case where 
pi and p2 are pure and finally for general pi and p2- 

We know already that the effects E„i of an optimal 
measurement are pure, i.e.. 



Em = a„jO„i)(0^ 



(18) 



Here, the states |0m) span the whole Hilbert space but 
they are in general not mutually orthogonal. Since we 
are dealing with qubits the projectors |0m)(0m| can be 
expressed by means of Bloch vectors r„i = (Om|f?|Om): 



|0,„)(0,„| = ^(l + f™ -ct) 



(19) 



where the components of the vector a are the Pauli op- 
erators <7x, <Jy and Uz- Thus the effects read 



Em = _p„i(l 



(20) 



with Pm '■= OLml'^- From the completeness of the effects 
Em = 1 we obtain the following constraints for the 
weights pm and the vectors Tm'- 



1 and ^ 







(21) 



The optimal POVM maximizes Q in Eq. (fT3|) . Express- 
ing the states pa := {2pi + p2)/3 and pb := (pi + p2)/2 
in terms of Bloch vectors ~ tr[(7pi] 



Pa = ^{t + fa-a), pb = ^{l + n-^), (22) 



Q reads 



1y- i^+rm-ra)^ 



+ rm ■ n 



(23) 



The quantity Q can be further simplified by expanding 
it in powers of the difference vector Ar = — rt , 



1 V- (1 + rm ■ in + Ar)) 



1 + ?Vn ■ f h 



\}ZPm ( 1 + rm ■ [n + 2A?^ + ^ 

\ m 



(rm ■ Arl' 



+ rm ■ rb 



( r^ ■ M 
+ rm-rb 



(24) 



where we have employed conditions (PT|) to obtain the last 
line from the second. Starting with this representation of 



Q, wc can further reduce the set of interesting POVMs 
for our task. 

Proposition 3. When determining the optimal 
POVM for qubits it suffices to consider POVMs with pure 
effects and maximally three outcomes where the vectors 
rm lie in the plane spanned by fa, fb and in the Bloch 
sphere. 

Proof. In fact the effects Em of any qubit-POVM V 
can be represented by means of Bloch vectors fm- 



E„ 



- Pm(l 



(25) 



where the p,„ , fm satisfy condition (j2ip . Positivity of 
Em implies < \\fm\\ < 1 with ||rm|j = 1 iff Em is pure. 

Now, let us assume without restriction of generality 
that the plane spanned by fa,fb and is the x-z plane. 
Starting from the vectors fm which describe the effects 
Em of the POVM V we consider their projection = 
fm\xz on the x-z plane. This results in a new POVM V 
with the effects 



E„ 



:p„i(l +qm ■ (?) 



(26) 



This POVM possesses in general not only pure effects, 
since the qm are in the generic case not normalized. How- 
ever, Eqs. [m are still valid, which implies that 
QiV) = QiV). To proceed, we can now decompose the 
non-pure effects Em into pure ones. This can for example 
be done by means of spectral decomposition: 



Em = Pt 



Ei 



Ei 



(27) 



Pm(l ± 



where the eigenvalues arc given by p^/^ 
j|(f,n!|). Please observe, that the Bloch vectors ±qm :~ 
±(fm/||(fm|| of the pure effects E'm'^ point into the x-z 
plane. Thus it suffices to consider POVMs with pure 
effects (but maybe many outcomes) in the x-z plane. Fi- 
nally, we can choose extremal POVMs. If the POVM 
with many pure effects in the x-z plane is not extremal, 
we can write it as a convex combination of extremal ones, 
which obviously have to have pure effects in the x-z plane 
again. The point is, that it is known that for extremal 
POVMs with four outcomes the projectors cannot lie in 
one plane [llj. This proves the claim. □ 
In the following we study the problem of optimal mea- 
surements beginning with the simplest case, the optimal 
von Neumann measurement for pure states pi and p2- 
Then we attempt to show that this von Neumann mea- 
surement leads to the least mean variance compared to 
all generalized qubit measurements. We increase the level 
of generality by studying the case of mixed states pi and 

P2- 



A. Optimal PVM for mixtures of pure qubit states 

Let us first consider pure states pi and p2 of a qubit, 
and find the optimal von Neumann measurement, i.e.. 



5 







fl 













r\ - r2 



aX a = kn with fc = 0, 1, 2, . . ., i.e., if the Bloch veetors 
r± of the projectors point into the same (opposite) di- 
rection as Ar. This is the case for measurements of non- 
degenerated observables O which commute with pi — p2- 
It remains to show, that this von Neumann measure- 
ment is also optimal among all possible POVM. This we 
will do in the next section, where we solve the problem 
for general pi,P2- 



^\^\ / B. Mixtures of mixed states 




FIG. 1: Schematic figure of the Bloch representation for the 
case where pi and p2 are pure. The Bloch vectors rh and 
Ar oc f 1 — r 2 are orthogonal, because of the circle of Thales. 
r denotes the direction of the considered measurement. In the 
calculation, it turns our that for the optimal measurement r 
is parallel to fi — r2 . See text for further details. 



the PVM {P+ = = 1 - \il>){tp\} for which the 

expected variance of A is least. The projectors {/+, P-} 
can be expressed by means of normalized Bloch vectors 
r_|_ = {ip\(j\ij:) = f* and r_ = — r, 



P± = i(l±fa). 



(28) 



This implies that the weights which appear in the param- 
eterization of the pure effects (pil)) are given hy p± = 1/2. 

Inserting the weights p± and the Bloch vectors f± into 
Q as given by Eq. yields 



1 / ^ {r-AfY 

4 I 1 - (f- ff,)2 



(29) 



Without loosing generality we assume that pi and p2 
are represented by the Bloch vectors fi = (0, 0, 1) 
and r2 = (sinfl, O,cos0), respectively. In this case 
ffc = (fl + f2)/2 = (sin^e; + (1 + cos6i)e;) /2 and 
Ar ~ {— sinOex + {I — cos9)ez) /6 are perpendicular, 
see Fig. [TJ Because of this orthogonality Q in (|^ re- 
duces to 




(Ar) cos^ a 



9 • 2 

rt sm a 



(30) 



where a is the angle between f and Ar, Ar = 11^^ I ^-nd 
''b = llfbll < 1- Q assumes its maximum value 



1 



1 + (Ar) 



(31) 



Which is the optimal measurement, if pi and p2 are 
mixed qubit states? In order to answer this question, 
we first determine the optimal projector-valued measure 
and then show that this leads also to a maximal value of 
Q for all POVMs. According to Eq. Q for a given 

PVM can be written as 



Q{a) 




(Ar) cos^(a) 
' cos^ /3 



1 



(32) 



where a is the angle between the Bloch vector f which 
represents one projector of the PVM and Ar and (3 is the 
angle between rand r),. Introducing the angle 7 between 
Ar and ff, and taking into account that (3 = a -f 7, one 
can directly verify that Q assumes its maximum at 



ao = i arccos 



cos 7 



^(r2-2)(l-cos(27)) + l. 



— 7 + fcTT with A; = 0,1, 



(33) 



where the positive sign is for < 7 < tt, and the negative 
sign has to be taken for — tt < 7 < (as it is in Fig. 2). If 
both states were pure (as in the previous section), we had 
7 = 7r/2 and hence ao = kn. For mixed states, however, 
this is in general not the case, hence the projectors of 
the optimal PVM-measurement do not correspond to the 
eigenvectors of pi — p2 (see Fig. 2). Only if both states 
pi and p2 have the same purity, (that is, tr(pi) = tr(p2) 
which implies ri = r2 4^ 7 = 0) the best measurement is 
again that of a non-trivial observable O which commutes 
with pi — p2- 

The deviation from the pure state case a = k-n is shown 
in Fig. 3. Physically, this deviation may be explained as 
follows. If p2 is relatively mixed compared to pi, then no 
observable does efficiently resolve the value of A for small 
A. This is simply due to the fact that the possible states 
with a small A are close to the maximally mixed state. 
So it is favorable to measure in the direction of pi which 
gives at least a high resolution for high values of A. 

Finally, we have to show that the optimal projection 
measurement is also the optimal POVM measurement. 
To do so, we prove that the Q-value for all POVMs is 
less or equal than its maximum value Q^™ for PVMs. 
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FIG. 2: Schematic figure of the Bloch representation for the 
case where pi is pure, but p2 not. Now, the optimal measure- 
ment is a measurement in the direction of r which is tilted 
towards ri. 




FIG. 3: Deviation of optimal measurement direction r from 
direction of fi — r2 specified by the angle ao as function of 
angle 7 between Af and ft for rb = 0.8. For states pi, p2 with 
the same entropy (7 = ±-k/2) the deviation is zero. If, say, 
state pi corresponds to less entropy r is tilted towards the 
Bloch vector of pi. For example compare the case depicted in 
Fig. [2I where n/2 < 7 < and thus < ao < 7r/2, indicating 
the named tilting. 



the optimal projective measurement, with ao given in 
Eq. (1221). That is, we have for aU POVMs Q < Qf„^f . 
Proof. We introduce the positive function 



a f{a) 



cos^(a) 



1 — cos^ (a + 7) ' 
then the claim can be expressed as follows: 

Q < ^ ^p„/(a„0(l - ncosPm) < /(ao) 

m 

/(ttm) 



(36) 



/(ao) 



(l-rbCOS/3™) < 1, (37) 



where := c^m + 7- Since /(a) < /(ao) for all a e 
[0, 2tt], it follows that 



/(«m) 

' /(ao) 



(l-r{,COS/3™) < ^pjn{l~n cos Pm) = 1 , 

m 

(38) 

which proves the claim. The right-hand side is equal to 
unity, because of the constraints on weights Pm and the 
Bloch vectors fm of the effects (|2T|). i.e., J^mP^n = 1 and 

LmPmCOS/Jm = 0. □ 



C. Decoherence as an example 

As an example, let us discuss the following problem. 
We consider a two-level atom under the influence of de- 
coherence. The task is, to estimate the decay rate from 
one copy of the state. The evolution of the density matrix 
may be described by the master equation 



d 1 
P= -z[H,p] +Cp, 



dt 



(39) 



where H = hioaz/'i- is the Hamiltonian of the atom, and 
the incoherent evolution is of the Lindblad form |12l | 



Cp 



(l-s)((T+o- p + pa+a -2(7 pa+) 



B 
"2 

- — s(ct^(T^P + pa^ — 2a^ pa^^ 
2C-B 



p - a^pa^). 



(40) 



For a general Qubit-POVM, (see Eq. Q is given by 



Q 



(Ar) cos^g^ \ 
-I- n cos(am + 7) / 



(34) 



whereas for the maximum Q for projection measurements 
is of the form: 



QZ" = - ( 1 



(Ar) cos^ ao 
1 - cos2(ao -|- 7) 



(35) 



with ao given by Eq. p3p . We can formulate: 
Proposition 4. For qubits POVMs can not improve 



Here, s = (0|poo|0) denotes the population of the exited 
state in thermal equilibrium, = ax i icTy , and B and 
C are the decay rates of the expectation values of and 



It can be straightforwardly shown that for the case 
C — B the time dependent density matrix is given by 



p{t) = Up{t)U^ 



-,-Bt 



p(0) + (l-e-^*) 



1-s 



(41) 



with U = e '^*/''. Note that this form is independent 
from the starting density matrix p(0). Now we can ask: 
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assuming we know s, t and p(0), how can we estimate the 
decay parameter B from a single copy of p{t)7 

Let us assume that wc know that at least B e 
[Q,B'^'^^]. Then the state under consideration is, in the 
rotating frame, of the form 



p(A) - Api + (1-A)p2, A = e-^*G[e 



-tB" 



.!]• (42) 



While B e [0,5'"^''] is equidistributed (i.e., p{B) = 
-^jj^max g^j]^ ^-j^ ^Yic parameter A is not. Its proba- 
bility density is given by 



9(A) = 



1 



(43) 



as can be checked by direct calculation. Employing the 
density q{X) we can calculate the mean variance using 
Eq. 



^p(TO)Var™(A) 

m 



l-e- 



2tB'° 



(tr[S„{(l 



xtB^ 



-)P1 + (^)P2}])' 



tY[E„i{xpi + (1 - X)p2}\ 



1-e 



-tB" 



(44) 



The remaining optimization problem is essentially the 
same as in the previous case. The only difference is that 
the weights of pi and p2 in the previous case [(2/3, 1/3) 
and (1/2,1/2)] are now replaced by more complicated 
expressions. Due to that, the discussion from above for 
the case where pi and p2 are pure, is no longer valid, 
since the orthogonality of ft, and Ar is not guaranteed 
anymore. However, the solution of the general case is di- 
rectly applicable, the only difference is the new definition 
of Ta and f;,. Hence, the calculations of Sec. IV can be 
applied, and finally Eq. ([55]) solves the problem of the 
estimation of decay parameter. 



and if O denotes a pure state we have Tr{0'^) = 1 which 
is equivalent to = 1. 

Now one can apply the results of the qubit case: First, 
one can argue as in Proposition 3 that the effects of the 
optimal POVM are linear combinations of Go,Gi and 
G2. Then, the previous section allows to compute the 
optimal ao and the corresponding vectors. Note that the 
normalizations 1 /d and \/ cP ~ d in Eq. (j45p are chosen 
in such a way, that all the formulae of the qubit case can 
be applied without modification. 

The drawback of this ansatz is that the obtained op- 
timal vectors are not guaranteed to correspond to valid 
POVMs. Namely, in contrast to qubits, for higher dimen- 
sional systems the condition ^ - rf = 1 does only imply 
that tr(O^) = 1, but not that O is positive. The con- 
ditions for positivity for Bloch vectors in higher dimen- 
sional systems are quite involved (isl . [l^ . However, if 
the resulting solution yields positive effects, the obtained 
solution is clearly the optimal one. We will discuss now 
important examples when this is the case. 

Let us first consider the case when pi and p2 have sup- 
port on the same two-dimensional subspace of Tid- This is 
for instance the case if pi = \ipi){ipi\ and p2 = \4>2){ip2\ 
are pure states. Then, by choosing an appropriate ba- 
sis of the two-dimensional subspace, we can assume that 
pi and p2 are real. Hence, Gi,G2 can be chosen as 
the Pauli matrices on the subspace, and the solution of 
the two qubit case can directly be applied. Moreover, 
if pi = |'0i)('(/'i| and p2 = \''p2){4'2\ the optimal mea- 
surement consists of a von Neumann measurement of an 
observable commuting with l^/;!) (■01 1 — \'ip2){4'2\- 

The other important case occurs if the first state is a 
pure state, pi = mixed with white noise p2 = l/d. 

Then we have to choose Gi = (1 — dji/') {ip\)/Vd'^ — 1 and 
G2 arbitrary. The optimal measurement is then a von 
Neumann-Liiders measurement with two effects: Pi = 
\^){'4>\ andP2 = l-K')(V'l- 



V. OPTIMAL MEASUREMENT FOR HIGHER 
DIMENSIONAL SYSTEMS 

Let us now discuss the case of higher dimensional sys- 
tems. That is, we consider a single copy of the state 
Pa ~ Api + (1 — A)p2 where the pt are density matrices 
acting on a d-dimcnsional complex Hilbert space Ti.d- 

A general solution for this case is quite complicated. 
However, for many important cases the problem can be 
solved as follows. Let us choose d^ operators Gi,i = 
0,...,d^ — 1 such that they form an orthonormal basis 
of the operator space. This means, they are hermitian 
and fulfill tr(GiGj) = 6ij. We can choose them in such a 
way that that pi as well as p2 can be written as linear 
combinations of Go = l/-\/d, Gi and G2. Any hermitian 
operator with trace one can be written as 

d^-i 

= -{t + Vd'^-d ^ r,G,;) (45) 



A. Verifying the production of entangled states 

The last example is similar to a task often occurring 
in experiments. Namely, one aims to produce an entan- 
gled state \'>p), however noise is added during preparation 
process. This situation may be modeled by writing the 
actual prepared state as 

piX) = m{^P\ + il-X)^t. (46) 

Now one would like to know whether the state is entan- 
gled or not. Our results deliver now a possible strategy 
for this decision: one may estimate A with the methods 
outlined above, and then apply separability criteria to the 
state p(Aest)- However, since only one copy is available, 
this does not allow to detect the entanglement unambigu- 
ously. It is interesting to see that this method is different 
from the standard method for many copies. Then, entan- 
glement witnesses allow the unambiguous detection, since 
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for them a negative mean value guarantees entanglement 
[isl . fioj Indeed, these are different observables: Taking 
the two-qubit case and \^) — a\Q\) + /?|10) the optimal 
witness is W = |00)(00| + |01)(10| + |10)(01| + |11)(11| 
[l7t which is different from the observable which leads to 
the best estimate of A. 



B. Optimal measurement for commuting states 

Another case where we can determine the optimal mea- 
surement to estimate the qudit state p\ in ^ is given 
when pi and p2 are commuting states, 

bi,P2]=0. (47) 

It turns out that measurements leading to the least mean 
variance are von Neumann measurements of observables 
O which commute with pi and p2, that is, they have 
the same eigenvectors as pi,p2- The proof can be ac- 
complished in two steps. First we prove that we can 
restrict the search for optimal measurements to the class 
of POVMs with effects commuting with pi and p2, then 
we argue that among these measurements the projection 
measurements yield the best estimation of px. 
Proposition 5. Let pi and p2 be commuting, [pi,p2] = 
and let the POVM V = {£',„} represent a measurement 
to estimate the state px = Xpi + (1 — A)p2- Then there is 
a projective measurement V = {F„i} with [F„i,pi] = 
for all m which satisfies Q{'P) = QiV). 

Proof. The commuting states pi and p2 are simultane- 
ously diagonalizable, 

Pi=J2^^\^)i^\ and p2='^s,\i){i\. (48) 

i i 

Thus Q{V), see Eq. (fT3)) . can be expressed as 

Q{V)^-J2 "'V/'^J' ''^ (49) 

with Crni (*|£'rn|i) Now, the positivc operators F„i := 
J2i^rni\i){i\ form a POVM, i.e., they satisfy the com- 



pleteness relation: 

= = l (50) 

It is easy to see that the commutative POVM V — {F„i} 
leads to the same value of Q as T'. Hence, we have only 
to consider POVMs which commute with pi and p2 to 
find a measurement which maximizes Q. In addition we 
already learned in section IIIII that an optimal POVM is 
distinguished by pure effects. Together with the commu- 
tativity of the effects it follows that the optimal measure- 
ment is a projective one. □. 
VI. CONCLUSION 

In conclusion, we have studied parameter estimation 
for quantum states, when only one copy of the state is 
available. For one qubit, we solved the problem by ex- 
plicitly constructing the measurement that minimizes the 
deviation between the true value of the parameter and 
the estimated one, using a Bayesian estimator. Further- 
more, we showed how the results from the qubit case can 
readily be used to solve this problem for important higher 
dimensional cases. 

Our work can be extended into several directions. 
First, one may look at higher dimensional systems, try- 
ing to find general solutions for this case. Here, it would 
be of great interest to find cases where, contrary to the 
qubit case, general POVMs allow for a better parameter 
estimation than von Neumann measurements. For prac- 
tical purposes, it may also be relevant to develop optimal 
measurement strategies for several, but a finite number 
of copies. 
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